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C3 ■ ABSTRACT 



£> | Galaxies are found to obey scaling relations between a number of observables. These 

relations follow different trends at the low- and the high-mass ends. The processes 

[^ ' driving the curvature of scaling relations remain uncertain. In this letter, we focus 

on the specific family of early-type galaxies, deriving the star formation histories of a 
complete sample of visually classified galaxies from SDSS-DR7 over the redshift range 

^ 0.01 < z < 0.025, covering a stellar mass interval from 10 9 to 3 x 10 M Q . Our 

sample features the characteristic "knee" in the surface brightness vs. mass distribution 
at M* ~ 3 x 10 10 Mq. We find a clear difference between the age and metallicity 
distributions of the stellar populations above and beyond this knee, which suggests a 
sudden transition from a constant, highly efficient mode of star formation in high-mass 
galaxies, gradually decreasing towards the low-mass end of the sample. At fixed mass, 
our early-type sample is more efficient in building up the stellar content at early times 
in comparison to the general population of galaxies, with half of the stars already in 
place by redshift z ~ 2 for all masses. The metallicity-age trend in low- mass galaxies 
is not compatible with infall of metal-poor gas, suggesting instead an outflow-driven 
relation. 

Subject headings: Galaxies: formation — Galaxies: evolution — Galaxies: stellar con- 
tent 



Introduction 



The conversion of gas into stars in galaxies and the role of feedback mechanisms has been 
one of the key aspects of extragalactic astrophysics over the past decades. Several baryonic pro- 
cesses regulate the star formation efficiency within dark matter haloes. Simulations of galaxy 
formation and semi-analytical models have shown that it is only within h aloes in a mass range 
around Mh a ,in ~ M? hnrk ~ 10 12 M^ where baryons can form stars efficiently (ICattaneo et al.1 12011 



Moster et al. 



2010 



Bouche et al. 



2010 



Guo et al. 



2011 



). Above this limit, 



gravitational shock heat- 



ing and AGN feedback suppress the gas accretion (JDekel fc Birnboimll2006l : iKeres et al.ll2009l . 120051 ; 
Birnboim fc Dekelll2003l : ICattaneo et alJl2009U201ll ). For galaxies within haloes with masses below 
~ 10 12 Mq, other processes ar e usually invoked to explain the star form ation suppression, for in- 



stance, reionization of the IGM (|Mamon et al.ll20ld : ICattaneo et al.ll201ll ). The energy liberated by 
supernova expl osions can eject the gas from haloes with circular velocity < 100 km s _1 , quenching 
star formation (JDekel &: Silkill986l ). 



The fraction of mass ac quired via mergers is also a function of stellar mass. For example, 
the semi-analytical models of iDe Lucia et al.l (|2006l ) show that the number of effective progenitors 
of galaxies with stellar masses < 10 11 Mq is less tha n two, while th i s num ber can be as large 
as five for galaxies with M* ~ 10 12 M Q . According to ICattaneo et al.l (|201ll ). the dependence of 
feedback and merger processes on stellar mass defines three galaxy formation regimes. Stellar 
mass ~ 10 11 Mq marks the transition between two dominant mechanisms: gas accretion (M* < 
10 11 M Q ) and gas-poor mergers (M* > 10 11 Mq). A third regime, set immediately below ~ 10 11 M Q , 
is characterized by the increasing contribution of a population built by gas-rich mergers. The 
contribution of these mass-depe ndent processes explain the well known dichotomy among early-type 
galaxi es (e.g. iKang et al.ll2007f ). supported by many observations, like e.g. a cha racteristic mass 



20041 ). 



scale (Kauffmann et al. 



20031 ) . and a well-defined mass-metallicity relation (see e.g. iTremonti et al. 



If the dichotomy originates from the mass-dependent role played by feedback, gas accretion, 
gas-rich and gas-poor mergers, we would expect to find signatures of these processes on the forma- 
tion history of galaxies with respect to stellar mass. In this letter, the star formation history of a 
sample of galaxies over a wide range of stellar mass (from 10 9 to 10 11,5 Mq) has been examined for 
the presence of those signatures. 



This letter is organized as follows: in Sect. [2j we describe the sample; in Sect. [31 we present a 
detailed study of the stellar populations using a spectral fitting code, which also is able to return 
the star formation history. Finally, we summarize and discuss our results in Sect. |U Throughout 
the paper, we adopt a cosmology with Hq = 70 km s _1 Mpc -1 , f2 m = 0.3 and Q\ = 0.7. 
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Fig. 1. — Upper panels: Scaling relations between stellar mass and surface brightness (left) 
or half-light radius (right). Both quantities are measured with respect to the Petrosian radius 
measured in the r band. Bottom panels: Mass- weighted average age (left) and metallicity (right), 
as obtained by the STARLIGHT spectral fitting code (see text for details). The grayscale corresponds 
to the density of points in each graph, with contour lines given at the 90, 70, 50, 25 and 5% levels 
with respect to the maximum value. The vertical dashed lines indicate M* ~ 3 x 10 10 M , which 
corresponds to the stellar mass of objects with .Mp e tro,T- = —20.2, i.e. at the knee in the scaling 
relations. Median values obtained within mass bins (see Table [1]) are shown in each panel as open 
dots. 



2. Sample description 



Our sample of early- type galaxies was retrieved from SDSS-DR7 (JAbazajian et al.ll2009l ). se- 
lecting galaxies in the redshift range between 0.01 and 0.025, brighter than m pc tro,r < 17.77, where 
r npetTo,v is the Petrosian magnitude in the r-band. This limit roughly corresponds to the magni- 
tude at which the SDSS spectroscopy is complete (jStrauss et al.ll2002l ). The redshift limits chosen 
provide a 95% complete sample between .Mpetro.r ~ —20 and ~ —17.46, where -M pc tro,r is the k- 
correct ed SDSS Petro s ian ab solute magnitude in r-band, obtained with the kcorrect code (version 
4_2) of Blanton et al. (200a), choosing as reference the median redshift of the sample (zq=0.021). 
( 2008 ) for details on the estimation of the completeness limits. 



See La Barbera et al. 



To select objects from the SDSS database, we define a flag mask, excluding those objects which 
are: i) blended, i. e., with more than one photometric pealo H) too bright (detections of > 200 ex); 
Hi) too large (r > 4' or a deblend with r > 1/2 frame); iv) saturated; v) located close to the edge of 
a frame or vi) in a region where the sky measurement failed, thus the photometry is compromised; 
objects which vii) are part of the extended wing of a bright star or viii) which may be an electronic 
ghost of a bright star were also excluded. In addition, only objects with no spectroscopic warning 
on (i.e. zWarning attribute set to zero) were selected. This selection returns 10187 objects. 

Early- type galaxies obey well-studied scaling relations. Several morphologic al classification 
indic ators have been proposed from the parameters of the SDSS pipeline (e.g. IStrateva et al 
200ll ). However, it is still unclear whether these indicators can be applied to low mass systems 



Previous work has shown that dw arf elliptical galaxies do n ot have the same surface brightness 
profile as their gian t counterparts (jGraham fc Guzman! I2003T). or similar star formation histories 
(jKoleva et alJ 12009 ). even featuring weak spiral structures (jLiskerl 120091 ) . Hence, we cannot use 
the standard SDSS attributes such as eClass or fracDev. Instead, visual inspection is the most 
reliable indicator of an ea rly- type morphology. We use the classification from the Galaxy Zoo 
project (JLintott et al.ll201ll ). selecting only galaxies classified as elliptical. We found 10138 (99.7%) 
galaxies from our sample in the Galaxy Zoo database. Among them, 1359 are classified as ellipticals 
(spirals and unreliable classifications are discarded). The completeness constraint - in -M pe tro,r ~ 
results in a final sample of 1328 objects. A further visual inspection was carried out by the authors 
to confirm the morphological classification. 



Stellar content 



Age, metalli city, stellar mass and veloc ity dispersion were derived using the spectral fitting 
code STARLIGHT (jCid Fernandes et al.ll2005l ). Before running the code, the observed spectra are 
corrected for foreground extinction and de-redshifted, and the models are degraded to match the 



Since we are using objects from the primary catalog, this selection is equivalent to exclude objects with two or more peaks, 
which were not dcblcndcd 
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Fig. 2. — Fraction of assembled stellar mass as a function of redshift. The panels show the fraction 
of mass normalized by the mass of galaxies at z ~ 0.1. Black triangles and solid (dashed) lines show 
median (mean) values for the present sa mple of early- type galaxies, w ith lcr confidence intervals 
being marked as dotted lines. Results by iPerez-Gonzalez et al.l (|2008l ) are plotted as (blue) solid 
circles with error bars, and dashed blue curves. 



wavelen gth-dependent resolution of the spectrum of each galaxy, as described in lLa Barbera et al 
(|2010bh . 



We u sed SSP models based on th e Medium resolution INT Li brary of Empiri c al Sp ectra 
(MILES, ISanchez-Blazquez et all 120061) . u sing the code presented in IVazdekis et al.l (J2010l ). us- 
ing version 9.1 ( Falcon-Barroso et alj|201ll ). They have a spectral resolution of ~ 2.5 A, almost 
constant with wavelength. A basis grid of 162 SSPs was selected, covering ages in the range of 
0.07 - 126 Gy r, and with [M/H] =1-1.71, -1.31, -0.71, -0.38, 0.00, +0.20). The model s use a 
Kroupal (j200ll ) Universal IMF with slope = 1.30, and isochrones by iGirardi et al.l (J200CJ). The 



stellar masses - computed within the fiber aperture - are extended to the full extent of the galaxy 
by computing the difference between fiber and model magnitudes in the z band. The stellar mass 
is then log(M*) = log(M*)' + 0.4 (m flbe r,z - m modchz ). 

We compare the results from STARLIGHT with different set-up parameters and different grids. 
We find no systematic trends, with differences typically within ±20%. A detailed study of how 
results are affected by changes on the set-up parameters of STARLIGHT and different SSP model 
assumptions will be given in a forthcoming paper (Trevisan et al., in pre p.). Variations in th e 



SSP optical colors due to different IMF shapes are very small (see e.g. IVazdekis et al.ll2010l ). 
Since STARLIGHT uses all the spectral information available, different IMFs are not expected to 
affect our results. Hence, i f a systematic change of the IMF is present (as suggested by, e.g. 



van Dokkum fc Conrov 



2011 



) from low- to high-mass early-type galaxies, the net result would be a 
change of the stellar mass that corresponds to the position of the knee, keeping the derived stellar 
population properties presented here unchanged. 



Results and discussion 



In this letter, we study the star formation histories of a sample of ~ 1300 visually-selected 
elliptical galaxies by means of a spectral fitting method. For each galaxy we determine the stellar 
mass, metallicity, age and star formation history (SFH). Fig. Q] shows both the photometric scal- 
ing relations (upper panels), along with the scaling of the derived average ages and metallicities, 
weighted in mass (bottom panels). The plots of surface brightness //p e tro and Petrosian radius vs. 
stellar mass provide similar information, considering that //p e tro is derived from Petrosian radius 
and luminosity, with the latter following the stellar mass of a galaxy. We bin the sample into five 
subsamples in stellar mass, indicated by blue dots in Fig. [1] Table Q] presents the median properties 
of the stellar populations for each bin. Regarding the age and metallicity scaling relation (bottom 
panels in Fig. [1]), a clear change in the slope of [M/H] vs. stellar mass is apparent at ~ 3 x 1O 1O M0, 
equivalent to an absolute magnitude of Mdcv,r ~ —20.2, below which the metallicity decreases 
linearly with mass. This corresponds approximately to the position of the knee seen in the photo- 
metric scale relations. The age distribution is more complex, with a homogeneous population of old 
galaxies for M* > 10 10 Mq, and an increased scatter towards younger ages with decreasing mass. 

It is not clear whether the dicho tomy in structural properties has the same o r igin as the 



stella r population properties (see e.g. iGraham fc Guzman! 120031 : iJanz fc Liskerj 120081 : 1 Cote et al 



2008). The fact that the knee in the photometric scaling relations has a counterpart in the stellar 
population properties might indicate the processes regulating the star formation also affect the 
structural properties of galaxies. For example, feedback mechanisms might affect galaxy sizes, 
since the gas is pushed out of galaxies by outflows and might be converted into stars at large radii, 
(e.g. iFan et al.lhoosl . hoiol : LDamjanov et alJboOfll ). 



4.1. Star formation history 



Besides the averaged age and metallicity, spectral fitting provides a wealth of additional infor- 
mation on the star formation histories of individual galaxies. For a given spectrum, a STARLIGHT 
run returns the contribution, as a percentage of mass, from each basis SSP. This distribution traces 
directly the star formation history. For each galaxy in the sample, we determine the "cumulative" 
mass fraction, i.e. the fraction of stars older than a given age, as a function of age. Then, we average 
the cumulative distributions over all galaxies within each mass bin. The age of the distribution at 
the 50th and 80th percentiles in stellar mass is presented in Table [TJ Galaxies with mass > 1O 1O M0 
form their stars early and over a very short period of time, with 80% of their stars being older than 
~ 11 Gyr. Galaxies in the low mass bins also have a very old stellar population. However, the time 
required to form 80% of the stellar mass is ~ 5 — 6 Gyr longer than that required by more massive 
galaxies. Hence, in early-type galaxies, downsizing should be interpreted as a more extended period 
of formation in low mass galaxies (instead of an overall later process of formation). 



We compare our results with iPerez-Gonzalez et al.l (|2008l ) . Their study is based on a sample 



of ~ 28000 objects of all morphological types at < z < 4 to constrain the evolution of the stellar 
mass content in galaxies as a function of redshift. We rebin our sample into three subsamples with 
masses ranging fr om log(M^) = 9.0 — 10.0, 10 -0 — 11.0, and 11.0 — 11.5, which correspond to the 



first three bins of iPerez-Gonzalez et al.l (J2008J) . Figure [2] shows the cumulative mass fractions as 
a function of redshift. In all three bins, early- type galaxi e s are formed in a much more efficient 
process, in contrast to the sample of lPerez-Gonzalez et al.l (|2008h . although the difference is more 
pronounced in the most massive bin. The large difference is due to the fact that their sample 
includes galaxies of all morphological types. In contrast, the analysis of the stellar populations 
of visually classified early- type galaxies at z < 1 , yield a short-lived and early process of star 
formation (JFerreras et al.l l2009bl ) , consistent with our findings. Furthermore, the present sample 
is susceptible to aperture effects, since all galaxies are observed through a fiber with fixed angular 
diameter. Table [1] reports the mean aperture for each mass bin. The aperture A is defined as the 
ratio between the radius of the SDSS fiber and the half-light Petrosian radius measured in the r 
band, A = -Rfibcr/-Rpctro50,r- The mean aperture varies from A ~ 0.5 in the first two bins to ~ 0.2 
in the more massive bin. Assuming that th e internal metallicity gr adient of early-type galaxies 
varie s from about —0.4 kpc -1 at high mass ( La Barbera et al.ll201ll ) to negligible at lowest mass 
(e.g. iKoleva et al.ll201ll ). the above variation of A would imply a change of ~ 0.16 in [M/H], i.e. 
smaller than that of ~ 0.4 seen for the range of masses in Table [TJ In addition, the typical internal 
scatter of stellar population properties is larger than systematics due to aperture effects. Since age 
gradients are generally small in (massive) ETGs, aperture effects are also likely not to drive the 
variation of galaxy age with stellar mass. This conclusion is further supported by our analysis of 



the w aveband dependence of the Fundamental Plane relation of bright ETGs (JLa Barbera et al 



2010al ). as we found that total (i.e. within an infinite aperture) metallicity and age do actually 
increase with stellar mass. 



Table 1: Stellar population properties as a function of stellar mass. 



log(M*) interval 


9.2-9.7 


9.7-10.2 


10.2- 10.7 


10.7- 11.2 


11.2-11.7 


Number of objects 


156 


293 


410 


360 


84 


L-weighted Age (Gyr) 


6.5 ±3.6 


8.8 ±3.3 


10.3 ±2.3 


10.2 ±1.6 


9.8 ±1.2 


M-weighted Age (Gyr) 


9.1 ±2.8 


10.3 ±2.4 


11.5 ±1.8 


11.6 ±1.4 


11.8 ±1.2 


50% of stars older than (Gyr) 


11.2 


12.4 


12.4 


12.4 


12.4 


80% of stars older than (Gyr) 


5.3 


8.2 


11.5 


11.6 


11.7 


L-weighted [M/H] 


-0.3 ±0.2 


-0.2 ±0.2 


0.0 ±0.1 


0.1 ±0.1 


0.1 ±0.1 


M-weighted [M/H] 


-0.2 ±0.2 


-0.1 ±0.1 


0.1 ±0.1 


0.2 ±0.1 


0.2 ±0.1 


-tt'fiber/ -tt-petro50,r 


0.51 ±0.19 


0.50 ±0.19 


0.42 ±0.15 


0.30 ±0.09 


0.21 ±0.05 



4.2. Constraints of feedback processes 

These results indicate that the processes regulating star formation in the low- and high-mass 
regimes leave different signatures on the SFH. Tab. Q] shows that galaxies within the three more 
massive bins have similar ages and metallicities, with 80% of their stellar mass formed at approxi- 
mately the same redshift. The ages and metallicities are almost constant from the third to the fifth 
mass bins. On the other hand, the low mass bins show a gradual decrease of age and metallicity 
with decreasing stellar mass. Galaxies in all bins have roughly half of their stars formed before 
redshift z ~ 2 — 3. For the most massive galaxies, the additional 30% of the stars is formed within 
~1 Gyr. Hence, the star formation in these galaxies after a redshift z ~ 1 can be considered "resid- 
ual". On the other hand, galaxies with mass ranging from log(M*) = 9.2 to 9.7 have 80% of their 
star formed only at z ~ 0.5. Our results indicate that massive objects form faste r and low mass 



sytem s have a more extended star formation history than suggested by models (e.g. iDe Lucia et al 



20061 ). 



The feedback mechanism commonly invoke d to explain the regulation of star formation in low 



mass systems is supernovae-driven winds (e.g. iLarsonl 1 19741 : iFerrara &: Tolstoy! [2000; iKeres et al 



xp 
u 



2009). The chemical enrichment of low-mas galaxies indicate that the process regulating the star 
formation in these galaxies should be strong enough to eject metals out of the galaxy. However, it is 
not able to quench the star formation completely. Alternatively, these systems may be continously 
being fuelled by infalling IGM gas at low metallicity. 

Figure [3] shows the mass- weighted [M/H] as a function of stellar population age. The metal 
enrichment is shown down to the age when approximately 80% of the stellar mass is formed_|. 
Galaxies in the three most massive bins form their stars very quickly, and 80% of their stars are 
older than 11 Gyr. For this reason, only two points are shown for these three bins. The figure 
shows a consistent trend of chemical enrichment, whereby the younger populations are more metal 
rich, a result that supports the idea that the mass-metallicity trend cannot be explained by the 
infall of metal-poor gas, requiring instead a feedback mechanism that preferentially removes metals 
from low-mass galaxies. Supernova feedback can account for these "metal-enhanced outflows". 
Since metals are formed in SN events, SN-driven winds are metal-enhanced with respect to the 
star-forming gas. Most of the metals mixed with t he hot gas are able to leave the galaxy, whereas 



only a small fraction of c ool ISM gas is lost (see e.g. iTremonti et al.ll2004l : lMac Low &; Ferrara 



1999 



Ferrara &; Tolstoy! l2000l ) . This scenario is compatible with the [M/H] vs. stellar mass relation 
shown in Fig. [TJ SN-driven feedback is expected to be dominant in low-mass systems, whereas 
feedback processes in massive galaxies leave a sharp truncation in the SFH, giving rise to different 
slopes above and below the characteristic mass. The position of the knee in the photometric scaling 
relations suggests that the interplay of these feedback mechanisms leaves an imprint on galaxy sizes 
and surface brightness. Aperture effect is not expected to affect these results, as galaxies in a given 



2 Notice that, for a given mass bin, the lowest Age in Fig. 3 does not correspond exactly to the value of t so % in 
Tab. 1, because of the discreteness of SSP ages in the STARLIGHT input basis (see Sect.[3]|. 
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Fig. 3. — Mass-weighted metallicity as a function of stellar population age, in bins of stellar mass. 
The metal enrichment is shown until 80% of the mass is assembled. For this reason, only two points 
at 11.2 and 12.6 Gyr are shown for the three most massive bins. This corresponds to the age of the 
two oldest SSPs used in the spectral fitting. The solid (dashed) lines indicate the median (mean) 
values. Error bars reflect mainly the discreteness of the SSP model grid; the metallicities available 
in the grid are spaced by ~ 0.4. 



mass range have similar rgber/^petro ratios. In addition, we verified that the trend is not driven 
by the age-metallicity degeneracy. We ran STARLIGHT on simulated spectra with no chemical- 



10 



enrichment and similar S/N and star-formation history as galaxies at lowest-mass, and we do not 
find any significant correlation between [M/H]Mass and Age. 

The SFH of high-mass galaxies indicate that their stars were already formed be fore z ~ 2. Sev- 
eral studies have shown th at high mass systems are assembled mainly via mergers ()De Lucia et al 
20061 : ICattaneo et al.ll201ll ). However, our results constrain the main epoch of growth via mergers 
to z > 2, unless the mergers proceed mainly via gas poor progenitors. Alternatively, these systems 
can be already in place at high redshift, as suggested by the weaker num ber density evo lution 



with redshift for the most massive galaxi es between z ~ 1.5 and 0. (see e.g. IConselice et al 



Ferreras et al 



Guo et al 



2009a 



2011 



Banerii et al. 



Dekel et al. 2009b 



201CJ). The cold mode of galaxy formation (JKeres et al 



2007; 



2005; 



a|) can account for the required high star formation efficiency at 



high redshift. However, it is still not clear how this process can result in the spheroidal morphologies 
we see in place at z ~ 1. 
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